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ABSTRACT 

The anomalous microwave emission detected in the Perseus molecular com- 
plex by Watson et al. has been observed at 11 GHz through dual orthogonal 
polarizations with the COSMOSOMAS experiment. Stokes U and Q maps were 
obtained at a resolution of ~ 0°.9 for a 30°x30° region including the Perseus 
molecular complex. A faint polarized emission has been measured; we find 
Q = -0.2% ± 1.0%, while U = -3.41} j% both at the 95% confidence level with 
a systematic uncertainty estimated to be lower than 1% determined from tests 
of the instrumental performance using unpolarized sources in our map as null 
hypothesis. The resulting total polarization level is 11 = 3.41} 9%. These are the 
first constraints on the polarization properties of an anomalous microwave emis- 
sion source. The low level of polarization seems to indicate that the particles 
responsible for this emission in the Perseus molecular complex are not signifi- 
cantly aligned in a common direction over the whole region, as a consequence of 
either a high structural symmetry in the emitting particle or a low-intensity mag- 
netic field. Our weak detection is fully consistent with predictions from electric 
dipole emission and resonance relaxation at this frequency. 
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1. Introduction 

Since the Cosmic Background Explorer (COBE) experiment first detected dust corre- 
lated microwave emission in its maps (Kogut et al. 1996), a significant effort has been made 
by the scientific community to understand the origin of this anomalous emission and charac- 
terize its properties. Further statistical evidence has been found in observations and analysis 
by Leitch et al. (1997), de Oliveira-Costa et al. (1997, 1998, 1999, 2002, 2004), Mukherjee 
et al. (2001, 2003), Lagache et al. (2003), Finkbeiner (2004), Finkbeiner et al. (2004) 
and Davies et al. (2005). Still, the underlying mechanism for this emission is a matter of 
discussion due to the lack of measurements about its properties, and different models have 
been proposed to explain its observed characteristics. 

The COSMOSOMAS experiment"*^ of the Institute de Astrofisica de Canarias (lAC) is 
now able to provide frequency and sky coverage needed to increase our understanding of the 
statistical and physical properties of the anomalous microwave emission ( Fernandez- Cerezo 
et al. 2006; S. Hildebrandt aZ. 2006, in preparation). Watson a/. (2005, hereafter W05) 
have recently presented a direct detection of rising-spectrum emission by COSMOSOMAS in 
the Perseus molecular cloud (R.A.=55°.4, decl.=31°.8; J2000) that is an order of magnitude 
higher than what can be explained with standard Galactic mechanisms of emission (i.e. , free- 
free, synchrotron, and thermal dust) and which cannot be explained with ultracompact H II 
regions or a gigahertz-peaked source. In this Letter, we present polarization measurements 
of this anomalous emission performed with the COSMOSOMAS experiment. It is worth 
stressing that the W05 detection refers to an extended region that appears diffuse even at 
the COSMOSOMAS angular resolution 1°) with a Gaussian FWHM fit over the emitting 
region of ~ 2°. Our polarization measurements are characterized by the same angular 
resolution. 



2. Anomalous microwave emission: interpretation and poleirization properties 

Comparisons between anomalous-emission detections and Ha maps have already ruled 
out the possibility that the statistical and direct observations are due to free-free emission 
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from ionized gas (see, e.g., Draine & Lazarian 1998a, 1998b). Ultracompact H II regions have 
also been invoked by McCuUough & Chen (2002) to explain the direct tentative detection 
obtained by Finkbeiner et al. (2002) as the superposition of a compact, optically thick 
and an extended, optically thin H II region, both being frec-frcc emitters. Bremsstrahlung 
emission is intrinsically unpolarized, although polarization may occur by means of Thomson 
scattering when photons are rescattered within an H II region. This may occur in optically 
thick regions, where the level of polarization can be at most 10% (Keating et al. 1998). 
Bennett et al. (2003a, 2003b) and Hinshaw et al. (2006) interpret the dust-correlated 
component as synchrotron emission with a flatter spectral index. However, the results of 
the Tenerife Experiment {e.g., de Oliveira-Costa et al. 2004) and the analysis of Lagache 
(2003) are not consistent with this interpretation. In any case, polarization measurements 
are essential to probe this hypothesis, because the synchrotron emission is expected to be 
highly polarized. 

Draine & Lazarian (1998a, 1998b) have proposed that the anomalous dust-correlated 
emission is due to electric dipolc emission from rapidly rotating small dust grains {i.e. "spinning 
dust" ) in the interstellar medium. They have delved into the details of the mechanisms of ex- 
citation and damping and concluded that their emission spectrum may fit well the observed 
signal and be responsible for this anomalous emission. Lazarian & Draine (2000) continued 
the study of spinning dust grains in the presence of weak magnetic fields and found that 
paramagnetic relaxation resonance, in a domain where classical paramagnetic relaxation is 
suppressed, may be efficient at producing an alignment of grains rotating faster than 1 GHz. 
This may result in the presence of an observable polarization as high as 5% at 10 GHz that 
rapidly decreases at higher frequencies, in the anomalous emission, depending on the phase 
in the sky of the magnetic field. Still, the level of alignment depends on factors such as the 
efficiency of spin-lattice relaxation which is uncertain for very small grains (< 10~^cm). 

Magnetic dipole emission from dust grains has also been found to be of considerable 
interest at frequencies below 100 GHz. Draine & Lazarian (1999) calculated the spectra 
of different grain candidates and found that ordinary paramagnetic grains exhibit emission 
spectra that are noticeably different from the observed anomalous emission. However, they 
proposed a model, adjusting the magnetic properties of the emitting material, that involves 
strong magnetic material whose emission may account, at least in part, for the observed 
anomalous emission. For this model, ferromagnetic relaxation may efficiently act to align 
dust grains and produce strongly frequency- and shape-dependent polarized emission that 
could be as high as 30% at 10 GHz. A key signature of polarized emission from a magnetic 
dipole, with respect to an electric dipole, is the variation of its direction with frequency. 

We note that both electric and magnetic dipole emissions are model dependent, making 
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it difficult to effectively predict their polarization properties. However, according to the 
model proposed by Draine & Lazarian (1999), in order to account for the observed spectrum, 
a hypothetical material "X4" characterized by strong magnetic grains has to be considered. 
Therefore, asymmetric grains will be aligned even by a low-intensity magnetic field like 
the Galactic field, resulting in both a higher polarization level and a higher probability for 
alignment to occur with respect to electric dipole emission. 

Iglesias-Groth (2005, 2006) studied the rotation rates of hydrogenated fullerenes and 
electric dipole emission in the interstellar medium and found that the smallest of these 
molecules could be the origin of the anomalous emission detected by W05 in the Perseus 
molecular complex and by Casassus et al. (2006) in the dark cloud LDN 1622. In addition, 
weak ferromagnetic properties exhibited by some of these molecules may lead to a consistent 
alignment and consequently to detectable polarization. 

3. COSMOSOMAS observations 

3.1. The instrument 

The COSMOSOMAS experiment is located at the Teide Observatory at 2400m above 
see level in Tenerife (Spain). It comprises two circular-scanning instruments and gener- 
ates two daily maps at 11 GHz (COSMOll instrument) and three maps at 13, 15, and 
16 GHz (C0SM015 instrument). The angular resolution is approximately 0°.9. The inter- 
ested reader may refer to Gallegos et al. (2001) for a general description of the C0SM015 
instrument and the adopted observational strategy. COSMOll has a similar instrumental 
setup to that of C0SM015. The optical configuration is identical (dimensions are scaled 
with wavelength) with a rotating primary flat mirror and a secondary parabolic mirror fo- 
cusing the radiation into a cryogenically cooled receiver. The detectors are low noise HEMT 
amphfiers cooled to 20 K. 

The main difference between the two experiments is that COSMOll is optimized for 
polarization measurements. An ortho-mode-transducer (OMT) is used to separate two or- 
thogonal polarizations after the radiation collector and before the HEMT, which are followed 
by further amplification stages. This configuration allows one to sample the sky and perform 

measurements of two Stokes parameters at a time: I and Q (or U, depending on the orienta- 
tion in the sky of the measured polarization planes). In order to completely characterize the 
linear polarization of the measured emission, the sampled polarization planes can be rotated 
by 45° by simply rotating the front end of the experiment, allowing one to switch the receiver 
sensitivity from Q to U (and vice versa). A detailed description of the instrument will be 
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presented in R. Hoyland et al. (2006, in preparation). 

3.2. Data Analysis 

COSMO 11 creates daily maps with full coverage in right ascension and ~ 20° coverage 
in declination. Raw data are stacked in scan collections in which the detected data, resulting 
from the adopted circular scan strategy, are saved as a function of time. Each scan is the 
result of an average over 30 cycles (equivalent to 30 s). Atmospheric emission is evident in 
the scans as a pseudo-sinusoidal modulation due to changing air masses within the circles 
covered on the sky by the instrument. This modulation has been removed by calculating 
a template and subtracting it from our scans. The template was calculated by averaging 
over 2 hr, masking bright sources and the Galactic plane, and iterating the procedure. The 
time over which the template was calculated was chosen by considering the final signal-to- 
noise-ratio as a function of time. This procedure allows the subtraction of slow day-night 
modulation, including 1/f noise as well as other fixed systematics such as optical pickup, 
which is in any event strongly reduced by under-illuminating the primary mirror. We note 
that this data processing differs from that adopted in W05, and by Fernandez- Cerezo et 
al. (2006), and is aimed at preserving the large scales of the emission {i.e. , to extend the 
window function to lower multipoles /). After the cleaning procedure, we proceed to the 
projection of the scans into the sky map, using bright sources to check for consistency in the 
pointing reconstruction. The final pixelization is 1/3° x 1/3°. After that we sum the daily 
maps into cumulated maps. The main calibration is performed with Cygnus A on a daily 
basis. Absolute flux is taken from the model by Baars et al. (1977). The instrumental beam 
is also calculated using the main calibrator for each day. The overall beam is recalculated 
using the cumulated map, and afterward we recahbrate all the maps with fixed FWHMs 
with Gaussian fits over the calibration source. 

In calculating the source emission, we take great care to account for contamination that 
might be due to the Sun or the Moon in the proximity of the target sources. We have 
adopted a conservative attitude, removing from the analysis those data in which the Sun 
is closer than ~30° to the regions of interest. Our conservative attitude originates from 
concern about possible (even very low) asymmetric pickup resulting in a polarized signal. A 
significant effort has been made to track the instrument's stability using radio sources and H 
11 region emission and polarization with time. Among other sources, wc have checked in the 
reconstructed map the coordinates of 4C 39.25 and 3C 345 outside the Galactic plane, and 
3C 84 and IC 405 in the Galactic anticenter, close to the Perseus region of interest. These 
sources have allowed us to test possible effects arising from the distance to the calibration 
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source in our maps, and to check that possible misahgnment effects in the optical system of 
the instrument are not affecting the local reconstruction of the microwave sky in the region 
of interest. Also, the relatively close H 11 region NGC 1499 (the California Nebula) has 
been used to test for possible effects and systematics arising from measurements of extended 
sources compared with pointlike sources. All of these sources have shown a good level of 
stability or slow variation compatible with intrinsic flux variations. NGC 1499 has also been 
carefully studied to check for spurious effects and as a null hypothesis test, since diffuse H 
II regions are expected to mainly emit free-free radiation, which is intrinsically unpolarized. 

Possible diffuse polarized synchrotron emission has been checked using the maps from 
WoUeben et al. (2006) ^. WoUeben et al. point out the presence of depolarized regions 
around H II regions. Fortunately, this does not have an effect on the polarized emission 
of the Perseus molecular cloud, which is known to be closer to us than NGC 1499. A 
detailed analysis of the WoUeben et al. polarization maps, and the extrapolation toward 
COSMOSOMAS frequencies, shows possible residual polarized signal lower than 1% in the 
NGC 1499 region, even assuming a spectral index as high as -2.7. We have thus finally 
re-calibrated our Perseus emission with NGC 1499. 

The coordinate convention adopted in this Letter to define the Q and U Stokes parame- 
ters is as follows. At every point on sky, the X-axis points north, and the Y-axis points east. 
In the first COSMO 11 configuration, we measure Iqo (intensity along the X-axis) and /goo 
(intensity along the Y-axis), while in the second orientation of the system we obtain /+450 
and I -450 . Our definition of the Stokes parameters is Q = Iqo — J900 and U = /450 — 1^450 . 
For each of the two configurations, we have two determinations of the intensity, which we 
denote Iq = Iqo + Iqqo and lu = /450 + /-45° • Finally, the total polarization degree is defined 
as n = y/Q^+lP/I. 

Spurious instrumental polarization induced by oblique rcficctions in our off-axis system 
has been calculated, in a regime of ordinary skin effect, using the method presented by 
Renbarger et al. (1998), and found to be totally neghgible compared with other systematics. 
The two channels have almost identical spectral response, resulting in a negligible systematic 
effect arising from the passband mismatch between observed calibrator and the Perseus region 
(i.e. , the effect is canceled out). Further instrumental effects have been tested for by rotating 
the frontend of the instrument by 90° and checking for consistency. 

Measurements of Q were taken over the period between 2004 March and 2005 May, while 
measurements of U are extracted from the measurements performed after this month. The 
time coverage is not uniform, as a result of contaminations, variation in the atmospheric 
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conditions, and instrumental failures, which have been reduced to a minimum thanks to 
continuous monitoring of the instrument's performance. The integration time for the U- 
measurements is smaller than that for Q, resulting in higher statistical errors. The residual 
systematics arc monitored by measurements of NGC 1499, using this H II region as a null 
test. This is done in two steps: by directly observing NGC 1499 on the Cygnus-calibratcd 
map, and then by considering the possibility of a faint polarization of the main calibrator 
which could result in an apparent polarization signal of the H II region. Tests on the WoUeben 
et al. (2006) maps allowed us to check for this possibihty and, after accounting for this, to 
constrain the systematics for both the Q- and U-directions at a lower level than the statistical 
uncertainties. In order to double-check this result we have carefully monitored 3C 84, whose 
emission is expected to be variable with time but whose polarization is expected to be far 
below 1% ^. We find a final polarization level lower than 1%, which sets our systematics 
level. 



4. Results and discussion 

The effective polarized emission from the Perseus region has been calculated through 
a maximum likelihood analysis using the measurements of the partial stacked maps. This 
likelihood was built using a multivariate Gaussian distribution. Measurements at different 
epochs were assumed to be uncorrelated, so the covariance matrix in this case is diagonal. 
We describe the data with a two-parameter model, one for the value (assumed to be constant 
in time) for the polarization of the Perseus region, and another for the calibration uncer- 
tainty. This latter parameter is marginalized over by following the analytical prescription 
given by Bridle et al. (2002), assuming a Gaussian prior for its value with 10% width (overall 
calibration uncertainty). Note that in the case of the fractional polarization, the calibration 
uncertainty does not enter. Once the likelihood curves are computed, the confidence limits 
are derived from the 0.025, 0.5, and 0.975 points of the cumulative distribution function. 
Thus, our parameter estimate is the median of the marginalized posterior probability distri- 
bution function, and the confidence interval encompasses 95% of the probability. In Figure 
1, we present the map of the Perseus region emission and the two difference maps of 90° 
polarization orientation describing the Q and U parameters. 

We find a small polarization in both orientations: for the 0°/90° orientation we obtain 
Q/I = -0.2% ± 1.0%, while for the 45°/ - 45° one we have U/I = -3.4+Jj% both at the 
95% confidence level, for a total polarization of 11 = 3.4l}'g%. The polarization in the 0°/90° 
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orientation is thus consistent with a null measurement within the systematic and statistical 
uncertainties. A polarization in U is possible, both from a statistical and a systematic point of 
view. From the considerations highlighted in the §1, although model-dependent uncertainties 
are widely present, the low level of polarization in the Perseus anomalous emission appears 
to be inconsistent with magnetic dipole emission of highly oriented particles in the presence 
of a magnetic field, either because of low field intensity in the region, or because of the 
highly symmetric nature of the emitting particles. Our results favor electric dipole emission 
as being responsible, the emission properties of which are even more model dependent but 
whose polarization is limited to ~ 5%. In particular, as described by Lazarian & Draine 
(2000), paramagnetic relaxation resonance may produce an observable polarized signal of the 
same order as that observed by COSMOSOMAS for an equivalent grain radius of < 10"'' cm. 
Our observations are limited both by the interpretation of the models and by our angular 
resolution. In fact, wc cannot constrain magnetically emitting grains characterized by weak 
magnetic properties or by spherical shape, and we cannot monitor structures characterized 
by angular dimensions smaller than our angular resolution in the observed region. This latter 
problem is, however, somewhat lessened by the fact that the measurement from W05 is over 
an extended region and so we expect anomalous-emission properties to characterize the entire 
region. Further polarization observations with higher sensitivity and angular resolution, as 
well as monitoring other frequencies of interest for the anomalous emission, will further help 
in understanding the mechanism responsible for this emission. 



5. Conclusions 

We have observed the Perseus molecular cloud in order to extract information about the 
polarization properties of the anomalous microwave emission observed by W05. A careful 
analysis and control of systematics lead us to the conclusion that this emission is character- 
ized by a low level of polarization: 11 = 3.4^} 9% at 95% confidence level, with systematic 
uncertainties limited to 1%. To infer the real emission mechanism for the anomalous emission, 
local physical properties of the observed regions should be further investigated. However, 
the weak detected polarization seems to support the electric dipole emission model, with res- 
onance relaxation, over the dipole magnetic emission hypothesis. Further information about 
the effective origin of this emission will be available with higher resolution, higher sensitiv- 
ity, multifrequency observations already planned at the Teide Observatory and, hopefully, 
planned for other instruments. 

We acknowledge the support of the lAC and the staff of Teide Observatory for the 
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Fig. 1. — Stokes I, Q, and U COSMOSOMAS 11 GHz maps of the region around the Perseus 
molecular complex. In the I map, the brightest source is the NGC 1499 H II region. The 
fainter source at R.A.=55°.4, decl.=31°.8 (J2000) is the observed (anomalously) emitting 
region. The source 3C 84 is also visible in the intensity map, with coordinates R.A.=49°.9 
and decl.=41°.5 (J2000). The color bar on the left refers to the I map, while that on the right 
refers to both Q and U maps. The maps have been smoothed with a 3 pixel boxcar. The 
inset in the U map refers to the faint detection achieved in the Perseus molecular complex. 



